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Abstract. This paper presents a nonconforming finite element approximation of tlie space 
of symmetric tensors with square integrable divergence, on tetrahedral meshes. Used for 
stress approximation together with the full space of piecewise linear vector fields for dis- 
placement, this gives a stable mixed finite element method which is shown to be linearly 
convergent for both the stress and displacement, and which is significantly simpler than 
any stable conforming mixed finite element method. The method may be viewed as the 
three-dimensional analogue of a previously developed element in two dimensions. As in that 
case, a variant of the method is proposed as well, in which the displacement approximation 
is reduced to piecewise rigid motions and the stress space is reduced accordingly, but the 
linear convergence is retained. 



1. Introduction 

Mixed finite element methods for elasticity simultaneously approximate the displacement 
vector field and the stress tensor field. Conforming methods based on the classical Hellinger- 
Reissner variational formulation require a finite element space for the stress tensor that 
is contained in H{div,Q;S), the space of symmetric n x n tensor fields which are square 
integrable with square integrable divergence. For a stable method, this stress space must be 
compatible with the finite element space used for the displacement, which is a subspace of 
the vector-valued function space. It has proven difficult to devise such pairs of spaces. 
While some stable pairs have been successfully constructed in both 2 and 3 dimensions, 
the resulting elements tend to be quite complicated, especially in 3 dimensions. For this 
reason, much attention has been paid to constructing elements which fulfill desired stability, 
consistency, and convergence conditions, but which relax the requirement that the stress 
space be contained in H{div, fl; S) in one of two ways: either by relaxing the interelement 
continuity requirements, which leads to nonconforming mixed finite elements, or by relaxing 
the symmetry requirement, which leads to mixed finite elements with weak symmetry. In 
this paper we construct a new nonconforming mixed finite element for elasticity in three 
dimensions based on tetrahedral meshes, analogous to a two-dimensional element defined in 



The space Hk of shape functions on a tetrahedral element K (which is defined in (3.1) 
below) is a subspace of the space V2{K; S), the space of symmetric tensors with components 
which are polynomials of degree at most 2. It contains Vi{K; S) and has dimension 42. The 
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degrees of freedom for a G S/< are the integral of a over K (this is six degrees of freedom, 
since a has six components), and the integral and linear moments of an on each face of K 
(nine degrees of freedom per face). For the displacements we simply take Vi{K,M?) as the 
shape functions and use only interior degrees of freedom so as not to impose any interelement 
degrees of freedom. See the element diagrams in Figure [1} We note that, since there are no 
degrees of freedom associated to vertices or edges, only to faces and the interior, our elements 
may be implemented through hybridization, which may simplify the implementation. See 
^ for the general idea, or for close to the present one. 




Figure 1. Degrees of freedom for the stress a (left) and displacement u 
(right). The arrows represent moments of an, which has three components, 
and so there are 9 degrees of freedom associated to each face. The interior de- 
grees of freedom are the integrals of a and m, which have 6 and 3 components, 
respectively. 

After some preliminaries in section |2| in section [3] we define the shape function space Y^k 
and prove unisolvence of the degrees of freedom. In section |4] we establish the stability, 
consistency, and convergence of the resulting mixed method. Finally in section |5] we describe 
a variant of the method which reduces the displacement space to the space of piecewise rigid 
motions and reduces the stress space accordingly. The results of this paper were announced 
in [13]. 

As mentioned, conforming mixed finite elements for elasticity tend to be quite complicated. 
The earliest elements, which worked only in two dimensions, used composite elements for 
stress [221 E]- Much more recently, elements using polynomial shape functions were de- 
veloped for simplicial meshes in two [10] and three dimensions [H H] and for rectangular 
meshes [H [H] . Heuristics given in [10] and [1] indicate that it is not possible to construct 
significantly simpler elements with polynomial shape functions and which preserve both the 
conformity and symmetry of the stress. Many authors have developed mixed elements with 
weak symmetry [^lElEZlESlESllMlElElElEnillSlIIS], which we will not pursue here. 
For nonconforming methods with strong symmetry, which is the subject of this paper, there 
have been several elements proposed for rectangular meshes [29l |30l EH [121 ES], but very 
little work on simplicial meshes. A two-dimensional nonconforming element of low degree 
was developed by two of the present authors in [TT]. As shape functions for stress it uses a 15 
dimensional subspace of the space of all quadratic symmetric tensors, while for the displace- 
ment it uses piecewise linear vector fields. A second element was also introduced in pi], for 
which the stress shape function space was reduced to dimension 12 and the displacement 
functions reduced to the piecewise rigid motions. In ^] Gopalakrishnan and Guzman devel- 
oped a family of simplicial elements, in both two and three dimensions. As shape functions 
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they used the space of all symmetric tensors of polynomial degree at most k + 1, paired with 
piecewise polynomial vector fields of dimension k, for k > 1. Thus, in two dimensions and 
in the lowest degree case, they use an 18 dimensional space of shape functions for stress, 
while in three dimensions, the space has dimension 60. Gopalakrishnan and Guzman also 
proposed a reduced variant of their space, in which the displacement space remains the full 
space of piecewise polynomials of degree k, but the dimension of the stress space is reduced 
to 15 in two dimensions and to 42 in three dimensions. However, their reduced spaces have 
a drawback, in that they are not uniquely defined, but for each edge of the triangulation 
require a choice of a favored endpoint of the edge. In particular, in two dimensions, the 
reduced space of [18] uses the same displacement space as the non-reduced space of [IT] , 
uses a stress space of the same dimension, and uses identical degrees of freedom, but the 
two spaces do not coincide (since the space of [TT] does not require a choice of favored edge 
endpoints). 

The element introduced here may be regarded as the three-dimensional analogue of the 
element in [TT]. Again, they have the same displacement space and the same degrees of 
freedom as the reduced three-dimensional elements of [T8], but the stress spaces do not 
coincide. Also, as in the two-dimensional case, our reduced space is of lower dimension than 
any that has been heretofore proposed. 



2. Preliminaries 

Let i7 C be a bounded domain. We denote by S the space of 3 x 3 symmetric matrices 
and by L^(f2;M^) and L^(f2;§) the space of square-integrable vector fields and symmetric 
matrix fields on Q, respectively. The space if (div, Q; S) consists of matrix fields r e L^(f2; S) 
with row-wise divergence, div r, in L^(fi; M.^). The Hellinger-Reissner variational formulation 
seeks {a,u) e if (div, fi; §) x L^{n;R^) such that 



(2.1) 



{Aa : T + div T ■ u)dx = 0, t e H{dw,Q;E>) 

/ diYa-vdx= / f-vdx, v e L'^{n;R''). 
Jn Jn 



Here a : r denotes the Frobenius inner products of matrices a and r, and A = A{x) :§—!■§ 
denotes the compliance tensor, a linear operator which is bounded and symmetric positive 
definite uniformly for x G fi. The solution u solves the Dirichlet problem for the Lame 
equations and so belongs to ii^(i7;]R^). If the domain Q is smooth and the compliance 
tensor A is smooth, then (cr, m) G ii^(f2;§) x H'^{Q;R^) and 

(2.2) ||a||i + ||M||2<c||/||o. 

with a constant c depending on Q and A. The same regularity holds if the domain is a 
convex polyhedron, at least in the isotropic homogeneous case. See [25] . 

We shall also use spaces of the form H^iVL^X) where X is a finite-dimensional vector 
space and is a nonnegative integer, the Sobolev space of functions i7 — > X for which all 
derivatives of order at most k are square integrable. The norm is denoted by || ■ H^^fc or || ■ 



To discretize (2.1), we choose finite-dimensional subspaces S/i C L^(f2;§) and Vh C 



L^(f2;]R^). Assuming that E/j consists of matrix fields which are piecewise polynomial with 
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respect to some mesh Th of fi, we define div/jT G L^(i7;]R^) by applying the (row-wise) di- 



vergence operator piecewise. A mixed finite element approximation of (2.1 ) is then obtained 
by seeking {ah,Uh) G S/j x \4 such that: 



/ {Aah : T + Uh- divh T)dx = 0, t E T^h 
Jn 

divh (Th - vdx = f ■ Vhdx, v e Vh- 
n Jn 



(2.3) 



If S/i C H{div,Q;S) this is a conforming method, otherwise, as for the elements developed 
below, it is nonconforming. We recall that a piecewise smooth matrix field r belongs to 
iJ(div) if and only if whenever two tetrahedra in Th meet in a common face, the jump IrnJ 
of the normal components rn across the face vanish. 



3. Definition of the new elements 

We define the finite element spaces S/i and Vh in the usual way, by specifying spaces of 
shape functions and degrees of freedom. The space Vh is simply the space of all piecewise 
linear vector fields with respect to the given tetrahedral mesh Th of Q (which we therefore 
assume is polyhedral). Thus the shape function space on an element K E Th is simply 
Vk = Vi{K] M^), the space of polynomial vector fields on K of degree at most 1. For degrees 
of freedom we choose the moments ft—)- J^^v ■ w dx with weights w G Vk- Since no degrees 
of freedom are associated with the proper subsimplices of K, no interelement continuity is 
imposed on Vh. The associated projection Ph : L^(fi;]R^) — t- Vh is the projection. 

To define the space we introduce some notation. If m is a unit vector, let : M'^ — > 
be the orthogonal projection onto the plane orthogonal to u. Then Qu is given by the 
symmetric matrix / — uu'. For a tetrahedron K, let Ak{K) denote the subsimplices of 
dimension k (vertices, edges, faces, and tetrahedra) of K. For an edge e G Ai{K) let Se be 
a unit vector parallel to e and let n+ and n~ be unit vectors each normal to one of the two 
faces of K which contain e. For a face / G A2{K) we denote by n/ its unit normal. We can 
then define the shape function space 

(3.1) EK = {aeV2{K;S) \ Qs^aQsJe e V,{e;S)ye e A,{K)}. 

For a G V2{K; §), Qs^crQs^ \e is a quadratic polynomial on e taking values in the 3-dimensional 
subspace Qs^^Qse of Thus the requirement that Qs^o'Qse \e belong to Vi represents 3 linear 
constraints on cr, and so dim T,k > 60 — 3 x 6 = 42. We shall now specify 42 degrees of 
freedom (linear functionals) and show unisolvence, i.e., that if all the degrees of freedom 
vanish for some a G Sx, then a vanishes. This will imply that dimS/^ < 42, and so the 
dimension is exactly 42. 

The degrees of freedom we take are: 

(3.2) aufvds, v e Vi{f; M^), / G A2(ii'), (36 degrees of freedom), 

(3.3) / adx, (6 degrees of freedom). 
Jk 

The following lemma will be used in the proof of unisolvence. 
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Lemma 3.1. Let fi and fj be the faces of K opposite two distinct vertices Vi and Vj and 
let e be their common edge, with endpoints Vk and vi. Given /3,7 G M, there exists a unique 
p E V2{K) satisfying the following four conditions (see Figure^: 

(1) p|eGPi(e), 

(2) p{vk) = 13, p{vi) = 7, 

(3) p\f^±L-Vr{f.),Pf^±L. V,{f,), 

(4) Jj^pdx = 0. 

Moreover pivi) = pivj) = 3(/3 + 7)/2. 




p{vi) = 7 

p\eeVi{e) 



Figure 2. The conditions of Lemma ^JJ 



Proof. For uniqueness we must show that if p G V2{K) satisfies (l)-(4) with /3 = 7 = 0, 
then p vanishes. Certainly, from (1) and (2), p vanishes on e, and then, using (3), p vanishes 
on fi and fj. Therefore p = cXiXj where Aj G Vi{K) is the barycentric coordinate function 
equal to on fi and 1 at fj, similarly for Aj, and c is a constant. Integrating this equation 
over K and invoking (4) we conclude that p does indeed vanish. 

To show the existence of p G V2{K), we simply exhibit its formula in terms of barycentric 
coordinates: 

p = ^Xl + {P + 7)A,A, + 7A,^ + + 7)(A? + A|) 

+ (-5/3 - 7)(A. + X,)Xk + (-/3 - 57)(A, + A,)A; + 3(/3 + 7)A,A,. 
That this function satisfies (l)-(4) follows from the elementary formula 

ai\ ■ ■ ■ ad+i\d\ 



X' 



ITI, a G 



IT {\a\+d)\ 

for the integral of a barycentric monomial over a simplex T of dimension d, which can be 
estabhshed by induction (see, e.g., |16]). □ 

We are now ready to prove the claimed unisolvence result. 



Theorem 3.2. The degrees of freedom given by (3.2) and (3.3) are unisolvent for the shape 



function space T.^ defined by (3.1 ).■ if the degrees of freedom all vanish for some a G T,^, 
then a = 0. 
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Proof. Let gi = grad Aj be the gradient of the ith barycentric coordinate function. Thus gi 
is an inward normal vector to the face fi with length 1/hi where hi is the distance from the 
ith vertex to /«. Note that any three of the gi form a basis for and that Y2i9i ~ 0- 

For a G T^k, define aij = dji = g^agj G V2{K). We shall show that if a G T,k and all the 
degrees of freedom vanish, then cr.y = on for all i ^ j. This is sufficient, since, fixing j 
and varying i, we conclude that agj = 0, and, then, since this holds for each j, that cr = 0. 

If e is an edge of the faces fi and fj of K, which may or may not coincide, then aij = 



g'iCrgj = g'iQsf^QsQj- Thus, from the definition (3.1) of the space Sj^, aij is linear on e. In 
particular, an is linear on each edge of fi. Thus p := crjj|j. is a quadratic polynomial on fi 
whose restriction to each edge of fi is linear. Therefore, on the boundary oi fi, p coincides 
with its linear interpolant, and, since a quadratic function on a triangle is determined by its 
boundary values, p is linear. Thus an is actually a linear polynomial on /j, and, in view of 



the degrees of freedom (3.2), we conclude that an vanishes on fi. 



For any pair k) of distinct indices (that is, 1 < /, A; < 4 and I ^ k), define 
(3.4) Pik = (^ijivk), (3ki = (Tij{vi), 

where i,j are the two indices unequal to / and k. Now aij G V2{K) is linear on the common 
edge e of fi and fj, and, because of the vanishing degrees of freedom of a, aij is orthogonal to 



Vi on fi and on fj and has integral on K. Therefore, by Lemma 3.1 applied with p = aij, 
it is sufficient to show that Pik and Pki both vanish in order to conclude that aij vanishes. 
In fact, we shall show that the 12 quantities Pik, corresponding to the 12 pairs of distinct 
indices, satisfy a nonsingular homogeneous system of 12 equations, and so vanish. 
The lemma also tells us that (Tijivj) = 3{(3ik + /3fc/)/2. Interchanging j and k gives 

3 

(7ik{vk) = ^{(3ij + (3ji). 



Also, by definition, 

(3.5) (3jk = aii{vk). 



Combining (3.4)-(3.5) gives 



3 

cTijivk) + (^ikivk) + (^iiivk) = ^iPij + l3ji) + il3ik + l3jk)- 

But aij + aik + an = —an, which vanishes on fi and so, in particular, at the vertex Vk. Thus 
we have established the equation 

(3.6) a(A, + Pji) + h{Pik + Pjk) = 0, 

where a = 3, 6 = 2. 

For each of the 12 pairs {i, k) of distinct indices, we let j and / be the remaining indices 



and consider the equation (3.6). In this way we obtain a system of 12 linear equations in 



12 unknowns. If we number the pairs of distinct indices lexographically, the matrix of the 
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system is: 
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Its determinant is 16(2a — 6)^6^(a + 6)^, as may be verified with a computer algebra package. 
In particular, when a = 3, 6 = 2, the system is nonsingular. Thus all the Pij vanish as 
claimed, and the proof is complete. □ 

Having established unisolvency, the assembled finite element space Hh is defined as the set 
of all matrix fields r such that t\k € for all K & Th and for which the degrees of freedom 
( 3.2[ ) have a common value when a face / is shared by two tetrahedra inTh- If r G S/^, then 
the jump \Tnf \ of rnj across such an interior face / need not vanish, but it is orthogonal to 
Vi{f; M.^). The normal component |?^jTnj] is, by the definition of the shape function space, 
linear on each edge of / so belongs to Vi{f), and thus 



(3.7) 



In'fTnfj 



on/. 



for any interior face of the triangulation. 



4. Error analysis 

In this section, we show that the pair of spaces S/,, Vh give a convergent finite element 
method. The argument follows the one given in |1]J for the two-dimensional case. As usual, 
we suppose that we are given a sequence of tetrahedral meshes Th indexed by a parameter h 
which decreases to zero and represents the maximum tetrahedron diameter. We assume that 
the sequence is shape regular (the ratio of the diameter of a tetrahedron to the diameter of 
its inscribed ball is bounded), and the constants c which appear in the estimates below may 
depend on this bound, but are otherwise independent of h. 

We start by observing that, by construction, 

(4.1) diYhJ^hCVh. 

The degrees of freedom determine an interpolation operator 11^ : H^{Q;S) — ?■ Hh by 
j^iUhT -r)n- vds = 0, vE Pi(/), / G A,i%), 

[ {UhT-T)dx = 0, KeTh, 
Jk 
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where Ak{Th) = UKeTh ^k{K). Since 



(div HfiT — div t) ■ V dx = — I (H/jT — r) : e{v) dx + (H/iT — r)n ■ v ds = 0, 
K Jk JdK 

for r G H^{K] §), t> G Vx, K E Th, have the commutativity property 
(4.2) diYhUhT = Ph div T, TeH\n;S). 



Since div maps iJ^(i7;§) onto L^(f2;]R^), (4.2) imphes that div^ maps E/^ onto Vh- An 
immediate consequence is that the finite element method system (2.3 ) is nonsingular. Indeed, 
if / = 0, then the choice of test functions t = ah and v = Uh imphes that ci/i = and then, 
choosing r with div/j t = Uh, we get Uh = 0. 

For the error analysis we also need the approximation and boundedness properties of the 
projections Ph and Uh- Obviously, for the projection, we have 

(4.3) \\v - Phv\\o < ch'^Wvlm, 0<m<2. 

Since Uh is defined element-by-element and preserves piecewise linear matrix fields, we may 
scale to a reference element of unit diameter using translation, rotation, and dilation, and 
use a compactness argument, to obtain 

(4.4) \\r-nhT\\o<ch^\\T\U m = l,2, 

where the constant c depends only on the shape regularity of the elements. See, e.g., [ID] for 
details. Taking m = 1 and using the triangle inequality establishes boundedness of Uh- 

(4.5) ||n;,r||o < c||r||i. 

The final ingredient we need for the convergence analysis is a bound on the consistency 
error arising from the nonconformity of the elements. Define 

(4.6) Ehiu, t)= [e{u) : r + divh T-u]dx, u G H^Q; M^), r G S?,, + H{dw, fi; §). 

Jn 

If r G //(div, Q; §), then Eh{u, r) = 0, by integration by parts. In general, 
Eh{u,T) = / TnK-uds= {rnfj-uds, 

where, again, [rn/J denotes the jump of rnj across the face /. Only the interior faces enter 
the sum, since u vanishes on dVL. Now rnj = Qnfi^nf) + {n'jrTnf)nf, so 



Eh{u,T) = ^ < lQnf{Tnf)j-uds+ lnjTnfj{n'fU)ds 
JlQnfirnf)j-uds, 



where the last equality follows from (3.7). 

We let Wh C \4 be the subspace of the displacement space Vh consisting of continuous 
functions which are zero on the boundary of Vt. In other words, Wh is the standard piecewise 
linear subspace of i/^(f2; M'^). For any r G S/^ the jumps, |Tn/], are orthogonal to Vi{f] M^), 
so Eh{w, r) = for any w G Wh- 
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Lemma 4.1. We may bound the consistency error 

(4.7) |E^(M,r)| < c/i(||r||o + /i||div;,r||o)||u||2, t G S^,, u E H^niW") H H^n^W"). 
Furthermore, for any p G H^{Q;S) 

(4.8) \Eh{u,Uhp)\ < ch'^\\p\\i\\u\\2, u e H\n;R^)nH'^{n;M.^). 

Proof. For any r G S/^ we have Eh{u,T) = E^iu — u{,t), where u{ G Wh is the piecewise 
hnear interpolant of u. Referring to the definition (4.6), we obtain 

\Eh{u,T)\ < c(|| div/,r||o||M - m^||o + lkllo||e(M - uj^)\\o < c/i(||r||o + h\\ div/j r||o) ||m||2, 

which is (4.7). For the second estimate we use that Eh{u, Hhp) = E{u — u{, Hhp) = Eh{u — 
u{, UhP — p), which imphes that 

Eh{u,Uhp) = y2 / div hO^hp - p) ■ {u - ui) dx + / {Uhp - p) : e{u - ui) dx. 
Ken •'^ •'^ 



Utihzing the estimate (4.4), the bound 

\Eh{u,Y\.hp)\ < c(|| divp||o||M - M^llo + Hn^p - p||o||e(M - m^)||o < c/i^||p||i||m||2 
is an immediate consequence. 



□ 



With these ingredients assembled, error bounds for the finite element method now follow 
in a straightforward fashion. 



Theorem 4.2. Let {a,u) be the solution of (2.1) and {ah,Uh) the solution of (2.3). Then 
(4.9) 



- CT/illo < c/i||m||2, 

diver — div/iO"/j||o < c/i™|| divcr||„ 
11^ ~ ^/illo ^ c/iIImIU. 



< m < 2, 



Furthemore, if problem (2.1) admits full elliptic regularity, such that the estimate (2.2) holds, 
then 



\u 



Proof. Subtracting the first equations of (2.1) and (2.3) and invoking the definition (4.6) of 
the consistency error, we get the error equation 



(4.10) 



/ [A{a - at) : T + {u - Uh) ■ div,, r] dx = Eh{u, r), r G S/^. 
Jn 



Comparing the second equations in (2.1) and (2.3), we obtain divhO'h = -P/idivcx, which 
immediately gives the claimed error estimate on diva. Using the commutativity (4.2), we 
find that divhiJ^hO' — Ch) = 0. Choosing r = UhC — ah in (4.10), we get 

/ A{cr - (Xh) ■■ {UhCT - ah) dx = Eh{u, Uha - ah), 
Jn 



which implies that 



W - ah\\\ < \W - ^hCr\\A + 2Eh{u,Ilha - ah) 



where ||t||^ := J At : rdx. Combining with (4.4) and (4.7) we conclude that 

llo" — O'hW < c/i(||(t||i + ||'U||2) < c/l||'u||2. 
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which is the desired error estimate for a. 

To get the error estimate for u, we choose p G i7^(n,S) such that div p = P^u — Uh and 



p||i < c\\PhU — Uh\\o- Then, in hght of the commutativity property (4.2) and the bound 
(|4.5 ), r := Hhp G S/^ satisfies div/jT = PhU — Uh and ||r||o < c\\PhU — Uh |o- Hence, using 



(4.1), (4.10), and (4.7), we get 



\\PhU - UhWl = / div/i r ■ {PhU - Uh) dx = divh t ■ {u - Uh) dx 
Jn Jn 

(4.11) = - / A{a - ah) ■■ rdx + Eh{u,T) < c{\\a - ah\\o + h\\u\\2)\\PhU -UhWo. 



This gives \\PhU — Uh\\o < ch\\u\\2, and then, by the triangle inequahty and (4.3), the error 
estimate for u. 

To estabhsh the final quadratic estimate for Hm — Uh||o in the case of full regularity, we use 
a duality argument. Let p = ^{w), where w G H^{Q; M.^) r\H'^{Q] M^) solves the problem 



divy4 ^ e{w) = PhU — Uh- It follows from (2.2) that 
(4.12) 



1 + \\w\\2 < c\\PhU - M;,||o. 



By introducing wj^ G Wh as the piecewise linear interpolant of w, we now obtain from (4.11 ) 
that 



\\PhU-Uh\\l = - / A{a - Gh) -.Hhpdx + Eh{u,Ilhp) 
Jn 

= - A{a - ah) ■■ {Uhp - p)dx + Eh{u, Uhp) - / {a - ah) : e{w - w{) dx, 
Jn Jn 

where the final equality follows since 

{a - ah) : e{wi) dx = - ^ / divh{a - ah) ■ w[dx + Eh{w[, a - ah) = 0. 



KdTh 



However, by utilizing (4.4), (4.8), the estimate for \\a — ah\\o given in (4.9), combined with 



the approximation property of the interpolant w^, we obtain from the representation of 
\\PhU — UhWi above that 

\\PhU - UhWl < c(/i^||p||i||m||2 + ||a - ah\\\\ t{w - w[)\\q) 

< c/i^||m||2(||p||i + ||w||2) < ch^\\u\\2\\PhU - UhWo, 



where we have used (4.12) to obtain the final inequality. This gives \\PhU — Uh\\o < ch \\u\\2- 



As above, the desired estimate for Hm— -u/Jlo now follows from (4.3 ) and the triangle inequality. 

□ 



5. The reduced element 

As for the two-dimensional element in [TT], there is a variant of the element using smaller 
spaces. Let 

T{K) = {v eVi{K; R^) \ v{x) = a + bxx, a,beR^}, 
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be the space of rigid motions on K. In the reduced method we take Vk '■= T{K) instead 
of Vk = Pi{K\M?) as the space of shape functions for displacement, so the dimension is 
reduced from 12 to 6. As shape functions for stress we take 

t^K = {t eT.k \ div/j r e T }, 



so dimEic = 36. As degrees of freedom for Y^k we take the face moments (3.2) but dispense 



with the interior degrees of freedom (3.3) 



Let us see how the unisolvence argument adapts to these elements, 
vanishing degrees of freedom, then divr G T{K), and for all v G T{K), 



T : e(v) dx + / rn v ds = 0, 



If r G with 



(div t)v dx = 

using the degrees of freedom and the fact that e{v) = 0. Thus divr = on and for all 
V G Vi{K;R^), 



(v) dx 



K 



K 



(div t)v dx + 

JdK 



rn V ds = 0. 



This shows that j^rdx = 0, so all degrees of freedom (3.3) vanish as well. Therefore the 
previous unisolvence result applies, and gives r = 0. 

A similar argument establishes the commutativity of the projection into (the analogue 
of ( |4.2 )), and the analogue of the inclusion (4.1) obviously holds. The space Tik still contains 
Vi{K] S) so the approximability (4.4) still holds, but the approximability of Vk is of one order 
lower, i.e., in (4.3) m can be at most 1. As a result, the error estimates given by (4.9) in 



Theorem 4.2 carry over, except that m is limited to 1 in the error estimate for div a. 
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